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I review the status of neutrino masses and mixings in the light of the solar and atmospheric neutrino 
data in the framework of two—, three- and four neutrino mixing. 



1 Indications for Neutrino Mass: 
Two Neutrino Analysis 

I ^] first review the present experimental sta- 
tus for solar and atmospheric neutrinos and 
the results of the different analysis in the 
framework of two-neutrino oscillations. 

1.1 Solar Neutrinos 
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The sun is a source of v' e s which are pro- 
duced in the different nuclear reactions tak- 
ing place in its interior. Along this talk I will 
use the v e fluxes from Bahcall-Pinsonneault 
calculations ff These neutrinos have been 
detected at the Earth by seven experiments 
which use different detection techniques iH: 
The chlorine experiment at Homestake, the 
water cerencov experiments Kamiokande and 
Super-Kamiokande (SK) and the radiochem- 
ical Gallex, GNO and Sage experiments. Due 
to the different energy threshold for the de- 
tection reactions, these experiments are sen- 
sitive to different parts of the solar neutrino 
spectrum. They all observe a deficit between 
30 and 60 % which seems to be energy de- 
pendent mainly due to the lower Chlorine 
rate. To the measurements of these six ex- 
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periments we have to add also the new results 
from SNO. They are however still not in the 
form of definite measured rates which could 
be included in this analysis. 

SK has also presented their results after 
1117 days of data taking ontl: 
- The recoil electron energy spectrum: SK 
has measured the dependence of the even 
rates on the recoil electron energy spectrum 
divided in 18 bins starting at 5.5 MeV. They 
have also reported the results of a lower en- 
ergy bin 5 MeV < E e <5.5 MeV, but its 
systematic errors are still under study and 
it is not included in their nor our analysis. 
The spectrum shows no clear distortion with 
X } lat = 13/(17do/). 

-The Zenith Angle Distribution (Day/Night 
Effect) which measures the effect of the Earth 
Matter in the neutrino propagation. SK finds 
few more events at night than during the day 
but the corresponding Day-Night asymmetry 
A D /N = -0.034 ± 0.022 ± 0.013 is only 1.3cr 
away from zero. 

In order to combine both the Day-Night 
information and the spectral data SK has also 
presented separately the measured recoil en- 
ergy spectrum during the day and during the 
night. This will be referred in the following 
as the day-night spectra data which contains 
2 x 18 data bins. 

The most generic and popular explana- 
tion of the solar neutrino anomaly is in terms 
of neutrino masses and mixing leading to os- 
cillations of v e into an active (y^ and/or u T ) 
or sterile neutrino, v s . In Fig. [I] I show the 
allowed two-neutrino oscillation regions ob- 
tained in our updated global analysis of the 
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Table 1. Best fit points and GOF for the allowed solutions for combinations of observables. 
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Figure 1. Presently allowed solar neutrino parame- 
ters for two-nautrino oscillations by the global analy- 
sis from Ref. Q. The plotted regions are 90%, 95 and 
99% CL. 



solar neutrino data i. We show the possible 
solutions in the full parameter space for oscil- 
lations including both MSW and vacuum, as 
well as quasi- vacuum oscillations (QVO) and 
matter effects for mixing angles in the second 
octant (the so called dark side). These re- 
sults have been obtained using the general ex- 
pression for the survival probability found by 
numerically solving the evolution equation in 
the Sun and the Earth matter valid in the full 
oscillation plane. In the case of active-active 
neutrino oscillations we find three allowed re- 
gions for the global fit: the SMA solution, the 
LMA and LOW-QVO solution. For sterile 
neutrinos only the SMA solution is allowed. 



For oscillations into an sterile neutrino there 
are differences partly due to the fact that now 
the survival probability depends both on the 
electron and neutron density in the Sun but 
mainly due to the lack of neutral current con- 
tribution to the water cerencov experiments. 
In Table @ I give the values of the parame- 
ters in these minima as well as the GOF cor- 
responding to each solution. There are some 
points concerning these results that I would 
like to stress: 

(a) Despite giving a worse fit to the observed 
total rates, once the day-night spectra data 
is included the LMA gives the best fit. This is 
mainly driven by the flatness of the spectrum 
and it was already the case with the last year 
data. 

(b) The GOF of the LOW solution has in- 
creased considerably as it describes the spec- 
trum data very well despite it gives a very 
bad fit to the global rates. Notice also that 
LOW and QVO regions are connected at the 
99 %CL and they extend into the second oc- 
tant so maximal mixing is allowed at 99 % 
CL for Am 2 in the LOW-QVO region. 

(c) As for the SMA the result from the cor- 
rect statistically combined analysis shown in 
Fig. 0and in Table |l| indicates that the SMA 
can describe the full data set with a proba- 
bility of 34% but it is now shifted to smaller 
mixing angles to account for the flatter spec- 
trum. 

(d) Similar statement holds for the SMA so- 
lution for sterile neutrinos. 
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Thus the conclusion is that from the statisti- 
cal point of view all solutions are acceptable 
since they all provide a reasonable GOF to 
the full data set. LMA and LOW-QVO solu- 
tions for oscillations into active neutrino seem 
slightly favoured over SMA solutions for os- 
cillations into active or sterile neutrinos but 
these last two are not ruled out. 



1.2 Atmospheric Neutrinos 

Atmospheric showers are initiated when pri- 
mary cosmic rays hit the Earth's atmosphere. 
Secondary mesons produced in this collision, 
mostly pions and kaons, decay and give rise 
to electron and muon neutrino and anti- 
neutrinos fluxes. Atmospheric neutrinos can 
be detected in underground detectors by di- 
rect observation of their charged current in- 
teraction inside the detector. These are the 
so called contained events. SK has divided 
their contained data sample into sub-GeV 
events with visible energy below 1.2 GeV and 
multi-GeV above such cutoff. On average, 
sub- GeV events arise from neutrinos of sev- 
eral hundreds of MeV while multi-GeV events 
are originated by neutrinos with energies of 
the order of several GeV. Higher energy muon 
neutrinos and antineutrinos can also be de- 
tected indirectly by observing the muons pro- 
duced in their charged current interactions in 
the vicinity of the detector. These are the so 
called upgoing muons. Should the muon stop 
inside the detector, it will be classified as a 
"stopping" muon, (which arises from neutri- 
nos of energies around ten GeV) while if the 
muon track crosses the full detector the event 
is classified as a "through-going" muon which 
is originated by neutrinos with energies of the 
order of hundred GeV. 

At present the atmospheric neutrino 
anomaly (ANA) can be summarized in three 
observations: 

- There has been a long-standing deficit of 
about 60 % between the predicted and ob- 
served Vfj,/v e ratio of the contained events i 



now strengthened by the high statistics sam- 
ple collected at the SK experiment 0. 

- The most important feature of the atmo- 
spheric neutrino data at SK is that it exhibits 
a zenith-angle-dependent deficit of muon neu- 
trinos which indicates that the deficit is larger 
for muon neutrinos coming from below the 
horizon which have traveled longer distances 
before reaching the detector. 

- The deficit for thrugoing muons is smaller 
that for stopping muons, i.e. the deficit de- 
creases as the neutrino energy grows. 

The most likely solution of the ANA in- 
volves neutrino oscillations. In principle we 
can invoke various neutrino oscillation chan- 
nels, involving the conversion of into ei- 
ther v e or v r (active-active transitions) or 
the oscillation of into a sterile neutrino 
v s (active-sterile transitions) 0. Oscillations 
into electron neutrinos are nowadays ruled 
out since they cannot describe the measured 
angular dependence of muon-like contained 
events 13. Moreover the most favoured range 
of masses and mixings for this channel have 
been excluded by the negative results from 
the CHOOZ reactor experiment i. 

In Fig. H I show the allowed neutrino os- 
cillation parameters obtained in our global fit 
of the full data set of atmospheric neutrino 
data on vertex contained events at 1MB, Nu- 
sex, Frejus, Soudan, Kamiokande □ and SK 
experiments □ as well as upward going muon 
data from SK, Macro and Baksan experi- 
ments in the different oscillation channels. 

The two panels corresponding to oscil- 
lations into sterile neutrinos in Fig. || differ 
in the sign of the Am 2 which was assumed 
in the analysis of the matter effects in the 
Earth for the — ► v s oscillations. Concern- 
ing the quality of the fits our results show 
that the best fit to the full sample is obtained 
for the — > v T channel although from the 
global analysis oscillations into sterile neutri- 
nos cannot be ruled out at any reasonable 
CL. Due to matter effects the distribution 
for upgoing muons in the case of — > v s 
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Figure 2. Allowed atmospheric oscillation parame- 
ters all for experiments, combined at 90 (shadowed 
area), 95 % and 99 % CL (thin solid-Line) for all pos- 
sible oscillation channels, from Ref. 0. The expected 
sensitivity for upcoming long-baseline experiments is 
also displayed. 



are flatter than for — > v T . Data show a 
somehow steeper angular dependence which 
can be better described by — ► v T . This 
leads to the better quality of the global fit 
in this channel. Pushing further this feature 
SK collaboration has presented an analysis of 
the angular dependence of the through-going 
muon data in combination with the up-down 
asymmetry of partially contained events and 
the neutral current enriched events which 
seems to disfavour the possibility — > v s 
at the 3-<7 level 

2 Three Neutrino Oscillations 

In the previous section I have discussed the 
evidences for neutrino masses and mixings as 
usually formulated in the two-neutrino oscil- 
lation scenario. Let us now fit all the different 
evidences in a common three-neutrino frame- 
work and see what is our present knowledge 
of the neutrino mixing and masses. Here I 
present a brief summary of such analysis per- 
formed in Ref. □ and I refer to that publica- 
tion for further details and references. 

The evolution equation for the three neu- 
trino flavours can be written as: 
.dv 



dt 



2E 



(1) 



where M v is the diagonal mass matrix for the 
three neutrinos and U is the unitary matrix 
relating the flavour and the mass basis. Hi nt 



is the Hamiltonian describing the neutrino in- 
teractions. In general U contains 3 mixing 
angles and 1 or 3 CP violating phases de- 
pending on whether the neutrinos are Dirac 
or Majorana. I will neglect the CP violat- 
ing phases as they are not accessible by the 
existing experiments. In this case the mix- 
ing matrix can be conveniently chosen in the 
form 

U = i?2 3 (M X fii 3 (01 3 ) X #12(012) , (2) 

where Rij is a rotation matrix in the plane 
ij. With this the parameter set relevant for 
the joint study of solar and atmospheric con- 
versions becomes five-dimensional: 

Amli . A ™L = Am§ 2 , 

C7q = 6*12 , Oatm = $23 ; Qreac = $13 , 

(3) 

where all mixing angles are assumed to lie in 
the full range from [0, n/2]. 

In general the transition probabilities 
will present an oscillatory behaviour with two 
oscillation lengths. However from the re- 
quired hierarchy in the splittings Am„ tm 3> 
Am| indicated by the solutions to the so- 
lar and atmospheric neutrino anomalies it fol- 
lows that: 

- For solar neutrinos the oscillations with 
the atmospheric oscillation length are aver- 
aged out and the survival probability takes 
the form: 



Am| 



p3y 

r ee.MSW 



713 



COS 



l3 r ee,MSW 



(4) 



where msw ^ s obtained with the modified 
sun density N e — > cos 2 9i^N e . So the analy- 
ses of solar data constrain three of the five in- 
dependent oscillation parameters: Am 2 ,!, #12 
and . 

- Conversely for atmospheric neutrinos, the 
solar wavelength is too long and the corre- 
sponding oscillating phase is negligible. As 
a consequence the atmospheric data analysis 
restricts Ato| 2 , 023 and 0i3, the latter be- 
ing the only parameter common to both so- 
lar and atmospheric neutrino oscillations and 
which may potentially allow for some mutual 
influence. 
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Figure 3. Allowed regions at 90% and 99% CL for 
the oscillation parameters Am2i and tan 2 (6*12) from 
the global analysis of the solar neutrino data in the 
framework of three— neutrino osciUations for different 
values of the angle 6*13 from Ref. □• 



Therefore solar and atmospheric neu- 
trino oscillations decouple in the limit #13 = 
0. In this case the values of allowed parame- 
ters can be obtained directly from the results 
of the analysis in terms of two-neutrino os- 
cillations presented in the first section. De- 
viations from the two-neutrino scenario are 
then determined by the size of the mixing 
#13. This angle is constrained by the CHOOZ 
reactor experiment which imposes an strong 
lower limit on the probability 



P, 



CHOOZ 



l-sin 2 (2Msin 2 (^§^) (5) 



> 0.91 at 90 % CL for Am| 2 > 10" 3 eV 2 . 
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The first question to answer is how the 
presence of this new angle affects the analy- 



fs of the solar and atmospheric neutrino data 
In Fig. U I show the allowed regions for the 
oscillation parameters ATO21 and tan 2 ($12) 
from our global analysis of the solar neu- 
trino data in the framework of three-neutrino 
oscillations for different values of the angle 
0i 3 . The allowed regions for a given CL 
are defined as the set of points satisfying 
the condition ^(Am^, tan 2 612, tan 2 013) — 
Xmin — ^X 2 (CL, 3 dof) where, for instance, 
Ax 2 (CL, 3 dof)=6.25, 7.83, and 11.36 for 
CL=90, 95, and 99 % respectively. The 
global minimum used in the construction of 
the regions lays in the LMA region and cor- 
responds to tan 2 0i3 = 0, this is, for the "de- 
coupled" scenario. Notice that the only dif- 
ference between the first panel in Fig. ^| and 
the active oscillations solution in Fig. || is due 
to the different numbers of dof used in the 
definition of the regions. The behaviour of 
the regions illustrate the "tension" between 
the data on the total event rates which favour 
smaller 0i3 values and the day-night spectra 
which allow larger values. It can also be un- 
derstood as the "tension" between the energy 
dependent and constant pieces of the electron 
survival probability in Eq. (||) . 

As seen in the figure the effect is small 
unless large values of 0i3 are involved. From 
Fig. H we find that as tan 2 0i3 increases all 
the allowed regions disappear, leading to an 
upper bound on tan 2 0i3 for any value of 
Am 2 ,!, independently of the values taken by 
the other parameters in the three-neutrino 
mixing matrix. The corresponding 90 and 99 
% CL bounds are tabulated in Table ||. 

As for the atmospheric neutrino data in 
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Fig. |^ I show the (tan 2 023, Am| 2 ) allowed re- 
gions, for different values of tan 2 813 from the 
global analysis of the atmospheric neutrino 
data. The upper-left panel, tan 2 #13 = 0, cor- 
responds to pure — > v T oscillations, and 
one can note the exact symmetry of the con- 
tour regions under the transformation #23 — > 
7r/4 — #23- This symmetry follows from the 
fact that in the pure — > v T channel matter 
effects cancel out and the oscillation probabil- 
ity depends on #23 only through the double- 
valued function sin 2 (2023). For non- vanishing 
values of #13 this symmetry breaks due to the 
three-neutrino mixing structure even if mat- 
ter effects are neglected We see that the anal- 
ysis of the full atmospheric neutrino data in 
the framework of three-neutrino oscillations 
clearly favours the — > v T oscillation hy- 
pothesis. As a matter of fact the best fit 
corresponds to a small value of #13 = 9°. 
With our sign assignment we find that for 
non-zero values of #13 the allowed regions be- 
come larger in the second octant of #23- No 
region of parameter space is allowed (even at 
99% C.L.) for tan 2 13 > 0.6. Larger values of 
tan 2 #13 would imply a too large contribution 
of — > v e and would spoil the description of 
the angular distribution of contained events. 
The mass difference relevant for the atmo- 
spheric analysis is restricted to lay in the in- 
terval: 1.25 x fCT 3 < Am 2 2 /eV 2 < 8 x fCT 3 
at 99 % CL. Thus it is within the range of 
sensitivity of the CHOOZ experiment and as 
a consequence the angle #13 is further con- 
strained when we include in the analysis the 
results from this reactor experiment. This 
is illustrated in Table|| where one sees that 
the limit on tan 2 O13 is strengthen when the 
CHOOZ data is combined with the atmo- 
spheric neutrino results. 

One can finally perform a global analysis 
in the five dimensional parameter space com- 
bining the full set of solar, atmospheric and 
reactor data. As an illustration of such anal- 
ysis I present in Table| the resulting bounds 
on 613. The final results from the joint so- 
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Figure 4. 90, 95 and 99% CL three-neutrino allowed 
regions in (tan 2 $23, Am| 2 ) for different tan 2 S13 val- 
ues, for the combination of global analysis of atmo- 
spheric neutrino data from Ref. The best-fit point 
is denoted as a star. 



lar, atmospheric, and reactor neutrino data 
analysis lead to the following allowed ranges 
of parameters at 99% CL 

1.1 x 1CT 3 < Am 2 2 /eV 2 < 7.3 x 1(T 3 (6) 
0.33 < tan 2 23 < 3.8 
tan 2 13 < 0.085 (if solar LMA) 
tan 2 0i3 < 0.135 (if solar SMA) 

In conclusion we see that from our statisti- 
cal analysis of the solar data it emerges that 
the status of the large mixing-type solutions 
has been further improved with respect to 
the previous SK data sample, due mainly 
to the substantially flatter recoil electron en- 
ergy spectrum. In contrast, there has been 
no fundamental change, other than further 
improvement due to statistics, on the status 
of the atmospheric data. For the latter the 
oscillation picture clearly favours large mix- 
ing, while for the solar case the preference 
is still not overwhelming. Both solar and 
atmospheric data favour small values of the 



additional #13 mixing and this behaviour is 
strengthened by the inclusion of the reactor 
limit. 

3 Four Neutrino Oscillations 

Together with the results from the solar and 
atmospheric neutrino experiments we have 
one more evidences pointing out towards the 
existence of neutrino masses and mixing: the 
LSND results. These three evidences can be 
accommodated in a single neutrino oscillation 
framework only if there are at least three dif- 
ferent scales of neutrino mass-squared differ- 
ences which requires the existence of a light 
sterile neutrino. Here I present a brief update 
of the analysis performed in Ref. E3 of so- 
lar neutrino data in such framework of four- 
neutrino mixing and I refer to that publica- 
tion for further as well as for the relevant ref- 
erences. 

In four-neutrino schemes the rotation U 
relating the flavor neutrino fields to the mass 
eigenstates fields is a 4x4 unitary mixing ma- 
trix, which contains, in general 6 mixing an- 
gles (I neglect here the CP phases). Existing 
bounds from negative searches for neutrino 
oscillations performed at collider as well as 
reactor experiments impose severe constrains 
on the possible mass hierarchies as well as 
mixing structures for the four-neutrino sce- 
nario. In particular they imply: 

(a) Four-neutrino schemes with two pairs of 
close masses separated by a gap of about 1 eV 
which gives the mass-squared difference re- 
sponsible for the oscillations observed in the 
LSND experiment, can accommodate better 
the results of all neutrino oscillation experi- 
ments. 

(b) In the study of solar and atmospheric 
neutrino oscillations only four mixing angles 
are relevant and the U matrix can be writ- 
ten as U — -R34 i?24 -R23 -R12 We choose so- 
lar neutrino oscillations to be generated by 
the mass-square difference between V2 and 
v\. With this choice the survival of solar 



zy e 's mainly depends on the mixing angle $12 
and it is independent of ^34. The mixing 
$23 and $24 determine the relative amount of 
transitions into sterile v a or active and v T 
only through the combination cos $23 cos $24 
(C23C24). We distinguish the following limit- 
ing cases: 

• c 23 c 24 = corresponding to the limit of 
pure two-generation v e — > v a transitions. 

• c 23 c 24 = 1 f° r which we have the limit of 
pure two-generation v e — > v s transitions. 

• If C23C24 7^ 1, solar VgS can transform in the 
linear combination v a of active and v T . 

In the general case of simultaneous v e — > 
v s and v e — > v a oscillations the corresponding 
probabilities are given by 

Pu t —*v, — c 23 c 24 (1 — P"e— »"e) 1 CO 
Pu^u a = (1 - C23C24) (1 - P Va ^u.) ■ (8) 

where Pu e —w e takes the standard two- 
neutrino oscillation for Am^ and #12 but 
computed for the modified matter potential 
A = Acc + c^c^Ajvc Thus the analy- 
sis of the solar neutrino data in the four- 
neutrino mixing schemes is equivalent to the 
two-neutrino analysis but taking into ac- 
count that the parameter space is now three- 
dimensional (Am^, tan 2 $12, cf 3 C2 4 ). 

I first present the results of the allowed 
regions in the three-parameter space for the 
global combination of observables. In Fig. || I 
show the sections of the three-dimensional al- 
lowed volume in the plane (Am^, tan 2 ($i2)) 
for different values of C23C24. The global min- 
imum used in the construction of the regions 
lies in the LMA region and for pure active 
oscillations value of cf^cfy = 0. As seen in 
Fig. |B| the SMA region is always a valid solu- 
tion for any value of c^c^. This is expected 
as in the two-neutrino oscillation picture this 
solution holds both for pure active-active and 
pure active-sterile oscillations. Notice, how- 
ever, that the statistical analysis is different: 
in the two-neutrino picture the pure active- 
active and active-sterile cases are analyzed 
separately, whereas in the four-neutrino pic- 
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ture they are taken into account simultane- 
ously in a consistent scheme. We see that in 
this "unified" framework, since the GOF of 
the SMA solution for pure sterile oscillations 
is worse than for SMA pure active oscillations 
(as discussed in the first section), the corre- 
sponding allowed region is smaller as they are 
now defined with respect to a common mini- 
mum. 

On the other hand, the LMA, LOW and 
QVO solutions disappear for increasing val- 
ues of the mixing c\^c\^. I list in Table || the 
maximum allowed values of c\^c\ 4 for which 
each of the solutions is allowed at a given CL. 
We see that at 95 % CL the LMA solution 
is allowed for maximal active-sterile mixing 
c 23 c 24 — 0.5 while at 99% CL all solutions 
are possible for this maximal mixing case. 
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Figure 5. Results of the global analysis for the al- 
lowed regions in Am| x and sin 2 #12 for the four- 
neutrino oscillations. The different panels represent 
the allowed regions at 99% (darker) and 90% CL 
(lighter). The best-fit point in the three parameter 
space is plotted as a star. 



Table 2. Maximum allowed value of c| 3 c| 4 for the 
different solutions to the solar neutrino problem 



CL 


SMA 


LMA 


LOW 


QVO 


90 


0.9 


0.44 


0.3 


forbidden 


95 


all 


0.53 


0.44 


0.28 


99 


all 


0.72 


0.77 


0.88 
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